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It is known that a keV scale sterile neutrino is a good warm dark matter candidate. We study how this 
possibility could be realized in the context of gauge extensions of the standard model. The naive expectation 
leads to large thermal overproduction of sterile neutrinos in this setup. However, we find that it is possible to 
use out-of-equilibrium decay of the other right-handed neutrinos of the model to dilute the present density of 
the keV sterile neutrinos and achieve the observed dark matter density. We present the universal requirements 
that should be satisfied by the gauge extensions of the standard model, containing right-handed neutrinos, to 
be viable models of warm dark matter, and provide a simple example in the context of the left-right symmetric 
model. 
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I. INTRODUCTION 

Dark matter (DM) is one of the experimentally observed in- 
dications of physics beyond the standard model (SM). A wide 
variety of astrophysical and cosmological observations con- 
firm that JIdm — 0.2 part of the total energy density of the 
Universe is composed of some form of nonbaryonic matter 
which interacts very weakly [1]. The most common parti- 
cle physics explanation comes in the form of weakly interact- 
ing massive particles, which are heavy and weakly interacting 
thermal relics, leading to cold dark matter. Another common 
candidate for cold dark matter is the axion, which is light, but 
due to a specific generation mechanism it has an extremely 
small temperature [2]. Hot dark matter has high velocities 
and a large free streaming length and it contradicts the experi- 
ment, because it prevents the formation of the observed small 
scale structures in the Universe. The intermediate situation, 
warm dark matter (WDM) is, however, less explored. It may 
even provide a solution to some of the problems of the DM 
simulations, reducing the number of Dwarf satellite galaxies, 
or smoothing the cusps in the DM halos. 

A natural candidate commonly considered for WDM is a 
light sterile neutrino [3]. 1 A simple realization is the ^MSM 
[6, 7], where only three singlet fermions, which have Ma- 
jorana masses and Dirac mixing with ordinary (active) neu- 
trinos, are added to the standard model. Then, the mass of 
one sterile neutrino can be chosen in the range of several keV 
and with very small mixing with the active neutrinos, it will 
provide a particle with the lifetime exceeding the age of the 
Universe, which can be the WDM candidate. The virtue and 
at the same time the problem of the model is, that the ster- 
ile neutrino with such a small mixing (the only interaction of 
this particle is via the Yukawa couplings) never enters into 
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the thermal equilibrium, and it can be produced only by some 
nonthermal mechanism. If this were not true and the neu- 
trino reached thermal equilibrium at some moment in the early 
universe, then without any additional mechanism, the ther- 
mal relics with mass of over about 90 eV would overdose 
the Universe. At the same time, one needs knowledge of the 
physics before the beginning of thermal evolution of the Uni- 
verse in order to calculate unambiguously the abundance of 
sterile neutrinos in the z^MSM (see Refs. [8-10]). 

The possibility analyzed in this article is opposite to the 
z^MSM. We assume, that there is some additional (gauge) 
scale between the electroweak and Planck scales, and that 
the sterile neutrinos are charged under these additional gauge 
transformations. 

It turns out that it is possible to reconcile the thermal over- 
production of the DM with the observations. To do this, the 
abundance of the sterile neutrino should be diluted after it 
drops out of the thermal equilibrium. This happens if some 
long-lived particle decays while being out of thermal equi- 
librium after the DM sterile neutrino freeze-out. This effec- 
tively reduces the amount of the DM sterile neutrino rela- 
tive to the overall energy balance of the Universe; see Fig. 1. 
It is also easy to find a candidate for this long-lived heavy 
particle — another (heavier) sterile neutrino in the model. We 
formulate the requirements on the properties of the DM sterile 
neutrino and the out-of-equilibrium decaying particle to make 
the model consistent with existing observations and bounds. 
This generic analysis, important for all possible models of this 
type, is made in the Sec. II. In the end of this section, all the 
requirements are summarized. 

There exist other ways to avoid the overproduction of the 
DM sterile neutrino in the analyzed class of the models, which 
we will only mention here. One possibility is realised if all the 
new gauge interactions are at the grand unified theory (GUT) 
scale, while the reheating after inflation leads to temperatures 
below the GUT scale. This situation is similar to the ^MSM, 
because the sterile neutrinos do not reach thermal equilibrium. 
Another possibility requires large (of the order of thousand) 
number of degrees of freedom at the moment of the sterile 
neutrino freeze-out, which does not seem natural. 

In Sees. Ill and IV, we analyze the possibility to realize 
these constraints in the simplest models. We then use other 
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under the gauge transformations of an extended model (ulti- 
mately, emerging in the breaking chain of some GUT model). 
Though for most of the statements in this article the precise 
details of this gauge interaction are not important, we will 
use a specific LR symmetric extension of the SM and stick 
to it to obtain definite numbers. This specific model (see 
e.g. Ref. [11] for a detailed review) with the gauge group 
5/7(3) x SU(2) L x SU(2) R x U{1) B -l appears as a sub- 
group of many GUT theories. 

In this model, we have the interaction with the gauge 
bosons of the form 



FIG. 1 : Schematic evolution of the light relic abundance in the Uni- 
verse. The dashed line is a thermal relic decoupled while being rela- 
tivistic (hot thermal relic), leading to the overclosure of the Universe. 
The blue decreasing line is the same hot thermal relic, but with the 
abundance diluted by rapid expansion of the Universe (entropy pro- 
duction), leading to correct DM abundance. The lowest (magenta) 
line depicts the evolution of the nonthermally produced particle with 
zero primordial abundance. 



sterile neutrinos to dilute the density of the DM sterile neu- 
trino. In Sec. Ill, we show that it is impossible for the same 
right-handed neutrinos to be involved in the DM abundance 
dilution and at the same time to give the masses to the active 
neutrinos via a type I seesaw like mechanism. The reason for 
this is that the mixing angles (or, equivalently, Yukawa cou- 
pling constants) are extremely small for the sterile neutrinos. 
This would lead to masses of the active neutrinos smaller, than 
the minimal ones, allowed by the neutrino oscillation obser- 
vations. 

In Sec. IV, we provide a working example, where the ac- 
tive neutrino masses are generated by a type II seesaw from 
the scalar sector of a left-right (LR) symmetric model, and 
sterile neutrinos have very small mixing angles with the ac- 
tive neutrino sector. 

The appendices are devoted to the calculation of the total 
decay width of the sterile neutrinos in the model (Appendix 
A), radiative decay width (Appendix B), and to the descrip- 
tion of the useful parametrization of the neutrino mass matrix 
(Appendix C). 



H. COSMOLOGICAL REQUIREMENTS AND 
CONSTRAINTS FROM EXPERIMENTS 

In this section, we introduce the generic framework we will 
work with and discuss the various constraints and bounds re- 
sulting from cosmological considerations and various experi- 
mental results. Note that these constraints are rather general 
and apply to most variations of the specified model. 



A. Assumptions and definitions 

In the following, we will assume the existence of right- 
handed (sterile) neutrinos Njr. These sterile neutrinos are not 
charged under the SM gauge group, but they could be charged 



- C CC = -j= ( W L laLl^aL + W& l a Rj„N aR ) 
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+ h.c. , (1) 

where Wl is the SM W boson, Wr is the corresponding right- 
handed boson from SU(2) r, and l a are the charged SM lep- 
tons. The neutrino mass matrix appears from the vacuum ex- 
pectation values of various Higgs bosons in the model. Up to 
the Sec. IV, we will not be interested in the details of this, and 
will just write the general mass matrix as 
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Note that a tilde over the neutrinos indicates that they are 
written in the flavor basis. In the following, we will also as- 
sume that the mass matrices obey in some sense the relations 
Mr > tud > Ml such that we can use seesaw-type formu- 
las. Thus, the rotation to the mass basis has the form 




(3) 



where U is the standard Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix and where Vr describes the mixing in the 
right-handed sector 

M L - mrjM^mJ = U* ■ diag(mi, m 2 , m 3 ) • <7 f , (4) 
Mr = • diag(Mi, M a , M 3 ) • V R , (5) 

with mi being the active neutrino masses and Mj the sterile 
neutrino masses. Note, that if Mr, = 0, then Eq. (4) is the 
usual seesaw formula. 

For the analysis of the sterile neutrino decay, when the os- 
cillations of the active neutrinos are not important, while the 
masses of the charged leptons are, it is helpful to make the 
described rotation only partially — without the PMNS rotation 
by the matrix U. Then, we get the mixing angles between the 
mass states of the sterile neutrinos and SM flavors 
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and also 



E 



(7) 



These squared mixing angles describe the overall strength of 
interaction (decay) of sterile neutrinos with the SM particles. 

Before moving on to the analysis of the cosmological prop- 
erties of sterile neutrinos, let us note an additional possible 
complication. Specifically, the Wl and Wr bosons in Eq. (1) 
may not coincide with the mass eigenstates, W\ and W% with 
masses M\y and M, respectively, but be slightly mixed 



Wl = cos C Wi + sin ( W 2 , 
Wr = - sin C Wi + cos C W 2 



(8) 



Normally this can be neglected, but it may give significant 
contribution to the radiative decay of the DM sterile neutrinos, 
analyzed in Sec. II F. 



B. Temperature of freezeout 

Let us now calculate the moment of decoupling of the neu- 
trinos Ni in the early universe. We will denote values corre- 
sponding to this moment by the subscript "f". As far as the 
DM sterile neutrino is relatively light and the freeze-out hap- 
pens while it is still relativistic, the calculation is analogous 
to those for the usual active neutrinos [12]. The only differ- 
ence is that the annihilation cross section is suppressed by the 
larger mass M of the right-handed gauge boson Wr, com- 
pared to the SM W boson mass Mw, 
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Here, ov P is the SM neutrino annihilation cross section, Gf = 
1.166 x 10 -5 GeV -2 is the Fermi constant and E is the en- 
ergy of the colliding neutrinos. Requiring the equality of the 
mean free path and the Hubble scale, we get for the freeze-out 
temperature 
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(10) 



where <j*f is the effective number of degrees of freedom im- 
mediately after freeze-out (at least 10.75 for SM content if 
freezeout happened below 100 MeV). 

We see that for the not very large scale M, the sterile neu- 
trino iVi decouples at a rather low temperature. Thus, it nor- 
mally is in thermal equilibrium at the early stages of the Uni- 
verse evolution, making it a thermal relic. This will be the 
possibility which we peruse in the current study. In this case, 
calculation of the present day density of the sterile neutrinos 
is insensitive to the history of the Universe before Tf. 

Note, however, that if the reheating temperature after infla- 
tion is lower than Eq. (10), the neutrinos never enter the ther- 
mal equilibrium. In this case, additional assumptions about 
the initial abundance of the sterile neutrinos are necessary to 



predict their current density, and the generation mechanism is 
very different form the analyzed here (see, e.g. Refs. [13-16]). 
Such a situation can be realized for a very low reheating tem- 
perature (see e.g. [17]), or naturally if the right-handed scale 
is the GUT scale, M ~ M GUT , leading to T f ~ M GUT , 
and the reheating after inflation reached slightly lower tem- 
peratures. Another way to implement this situation is perused 
in the z/MSM [6, 7], where no new physics is present up to 
Planck scale, leading to N\ never entering the thermal equi- 
librium. 



C. Abundance of Ni at present time 

The number to entropy density ratio of the sterile neutrino 
(two fermionic degrees of freedom) after freeze-out is given 



by 
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While the Universe expands slowly with all the processes ap- 
proaching thermal equilibrium, both the number density and 
entropy density scale are inversely proportional to the volume 
of the Universe, and this ratio remains constant. If nonequi- 
librium processes happen during expansion (for example an 
intermediate matter dominated stage caused by out of equilib- 
rium decay of a heavy species), additional entropy release is 
possible, which we will take into account by the factor S: 
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Let us calculate the contribution of N\ to the present en- 
ergy density. Rescaling the number to entropy density ratio 
at present moment by this factor, as compared to the freeze- 
out moment, we get for the sterile neutrino contribution to the 
energy density of the Universe fl jy 
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where f^DM = 0.105/i~ 2 is the DM density, s = 
2889.2 cm~ 3 is the present day entropy density, and p c = 
1.05368 x 10~ 5 h 2 GeV cm~ 3 is the critical density of the 
Universe. The observational requirement is Ojv/^dm < 1 
with equality being the nicest choice (all DM is made out of 
Ni) and inequality opting for multispecies DM. 

Let us analyze Eq. (13) further. Without entropy release 
(S = 1), the Universe is overdosed, unless the neutrino is 
very light, which corresponds to the hot dark matter case, ex- 
cluded by the structure formation in the Universe. Models 
with the number of degrees of freedom at freeze-out g„i of 
order 1000 seem rather unnatural and will not be considered. 
The only opportunity is thus the entropy release after freeze- 
out of Ni, 
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Having this entropy release after N± decoupling will lead to 
the observed DM abundance today. In the following, we will 
analyze possibilities of generation of this large amount in the 
model. 2 



D. Mass bounds 

The mass of the DM particle can not be too light, or the 
observed structure in the Universe would have been erased by 
a too hot DM. The simplest and most robust bound can be 
obtained from the phase space density arguments. The phase 
space density of a collisionless DM can only become smaller 
during the evolution of the Universe, as an effect of coarse- 
graining. Comparison of the primordial phase space density, 
which is calculated using the initial DM particle distribution 
function and the maximal modern one, derived from the ob- 
servation of the Dwarf spheroidal galaxies [4, 18] gives the 
lower bound 



Mi > 1-2 keV 
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Another important bound is the Lyman-a (Ly-a) bound 
[19, 20]. This bound constrains the velocity distribution of 
the DM particles from the effect of their free streaming on the 
formation of the structure on the scales, probed by the Ly-a 
forest. It should be noted that to convert this constraint into a 
bound for the mass of the DM particle, one needs to take into 
account the initial velocity distribution of the particles. In our 
case it takes the form of a usual thermal distribution, but with 
the temperature lowered by the dilution factor S^ 1 ^ 3 . This 
corresponds to the thermal relic case in Ref. [19], and not to 
the case of the nonresonantly produced sterile neutrinos, de- 
noted ttinrp in Ref. [19]. Thus, the result of Ref. [19] should 
be rescaled as 
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where T is the present temperature of the DM neutrino diluted 
with the entropy factor (14), T v is the temperature of the usual 
relic neutrinos, tonrp = 8 keV [19], and the ratio of the tem- 
peratures (T/Tv) 3 = f2 DM /i 2 (94 eV/Mi) is obtained from 
the requirement of the observed Odm- 



E. Generation of entropy 

The entropy (14) can be generated by some heavy long- 
lived particle which goes out of the thermal equilibrium at 
some moment after DM neutrino freeze-out tf, and it decays 
after becoming nonrelativistic and dominating the Universe 



2 The exact value of the required entropy release S may be slightly different 
if, for example, some amount of DM sterile neutrino was generated non- 
thermally after the freezeout. In the examples analyzed in the paper, this 
effect is negligible. 



expansion. The obvious candidates for such particles are the 
two remaining heavier neutrinos (though other candidates are 
possible and can be analyzed in a similar way). Let us assume 
for simplicity that only one of these two neutrinos is respon- 
sible for entropy generation, and we denote it by N%. Then, 
according to Refs. [12, 21], the entropy release is 
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where M 2 is the mass of N2, r = un 2 /s is the initial abun- 
dance of the N2 particles after decoupling (or, probably more 
precise, before they start to drive the matter dominated inter- 
mediate stage of the Universe expansion), and is the prop- 
erly averaged effective number of degrees of freedom during 
the N 2 decay. The ratio r is maximal when the particle decou- 
ples when it is still relativistic, and is equal to 
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where gpr = 2 is the number of degrees of freedom for N2, 
and g* is taken at the N2 freeze-out. 

If the entropy generation is large, we can neglect the 1 in 
Eq. (17) and get 
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By combining Eqs. (14) and (19), we obtain 
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Note that for our case, the freeze-out temperatures of the DM 
sterile neutrino and of the entropy generating one coincide, so 
g* = g*i- If Eq. (20) is satisfied, then we have proper DM 
abundance in the present Universe. 

However, Eq. (20) is not the only requirement for the life- 
time of the heavier sterile neutrino for the realistic model. En- 
tropy generation should finish before the big bang nucleosyn- 
thesis (BBN), i.e. N2 should decay before it. According to 
Refs. [22-24], the temperature after the decay of the sterile 
neutrino N2 should be greater than 0.7 4 MeV in order not 
to spoil BBN predictions. This temperature is approximately 
equal to (see Ref. [21]) 
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leading to the bound on the N2 lifetime which should be 
shorter than approximately 0. 1 -f- 2 s. The neutrino with such a 
lifetime can produce enough entropy, satisfying Eq. (20) only 
if it is sufficiently heavy, 
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Finally, as far as we were assuming that the sterile neutrino N2 
decoupled while still relativistic (otherwise the entropy gen- 
eration is much less efficient), we should require Tf > M 2 . 
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FIG. 2: Unitary-gauge diagrams contributing to the radiative neu- 
trino decay with charged leptons propagating in the loop. 



This, using Eq. (10), is translated into a bound for the scale of 
the right-handed bosons, 
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Thus, on the one hand the sufficient entropy generation re- 
quires a long-lived neutrino, but on the other hand, the require- 
ment of the successful BBN limits its lifetime from above, 
leading to the lower bounds on its mass and on the mass scale 
of the additional gauge interactions. 



F. x-ray observations 

A sterile neutrino in the considered class of the models is 
unstable, so it provides a decaying DM. Through its mixing, 
it decays via the neutral current into three active neutrinos. 
To lead to a successful DM scenario, the lifetime of the un- 
stable neutrino N\ should be greater than the age of the Uni- 
verse t u ~ 10 17 scc, which constrains its total decay width. 
However, one obtains significantly stronger restrictions result- 
ing from a subdominant decay channel — the radiative decay 
Ni —> vj, induced at the one loop level (Fig. 2). This process 
produces a narrow line in the x-ray spectrum of astrophysi- 
cal objects [3, 25]. In the context of ^MSM, the only source 
of this decay is via the active-sterile neutrino mixing 6\, and 
the recent x-ray observations bound it from above. A very 
rough bound, which will be enough for our purposes, is given 
inRef. [13] 3 
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This bound corresponds to the following upper bound on the 
radiative decay width 
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3 We must note that careful analysis gives a stronger (in some regions of 
masses by an order of magnitude) bound. See the detailed discussion in 
Sec. 5.1.2 of [13] and [26-32]. For our purposes, this approximate (weak) 
bound is sufficient. 



We also note that there are bounds resulting from supernova 
cooling. They are also much weaker than the diffuse x-ray 
background limits (24) for all possible neutrino masses Mi 
[3]. 

In the LR symmetric model, the x-ray bound (25) leads not 
only to the bound on the mixing angle (24), but also bounds 
the properties of the bosonic sector of the theory. The reason 
is the possible mixing of the right Wr gauge bosons with the 
SM Wl ones. Without mixing, the contribution of the Wr 
bosons to the process N — > ju is additionally suppressed by 
the ratio of the left and right gauge boson masses (Mw /M) 4 , 
and can be safely neglected. With the mixing, however, the 
chiral structure of the diagram changes, and the contribution 
can be enhanced by the factor (m;/Mi) 2 , where m; is the 
mass of the charged lepton running in the loop. 

We calculate the total decay width for Ni —> ^v, summed 
over the active neutrino flavors, following Refs. [33, 34] (for 
details, see Appendix B). Supposing from the very beginning 
that the right-handed scale is much larger than the left one, 
M 2> Mw, neglecting the active neutrino masses and assum- 
ing small gauge boson mixing, we get 
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is the mass 
The second 



Here, a is the fine-structure constant, and m; 
of the charged lepton propagating in the loop, 
term in the amplitude is proportional to the mass of the sterile 
neutrino Mi, while the first term to the mass of the charged 
intermediate lepton m; a . This can be understood from the 
following consideration. Because the photon has spin one, 
there must be a chirality flip on the fermionic line. If, in flavor 
basis, there is a Wr-Wl mixing, we have the chirality flip 
(mass insertion) on the internal line of the charged fermion, 
which produces a term proportional to mi a . Otherwise, the 
chirality flip happens on one of the outer lines of the diagram, 
with a term proportional to the Majorana mass of the incoming 
sterile neutrino. 

If the gauge boson mixing is absent, £ = 0, only the second 
term contributes and we obtain the usual result 
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If the mixing ( is present, then, barring the unlikely cancel- 
lation between the two terms in Eq. (26), we can constrain ( 
using the x-ray bound (25) 
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Thus, the mixing angle of the W bosons must be vanishingly 
small. 

We would like to note that in a LR symmetric model with 
the Higgs sector as described in Sec. IV, the Wl-Wr mixing 
angle £ is given by (see Ref. [11], and references therein) 
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where K\ and k-i are bidoublet vacuum expectation values 
(VEVs) and vji is the VEV of the right Higgs triplet. There- 
fore Eq. (28) restricts the ratio KiK^/v^. 



Note that this is the only requirement which changes 
in the case of entropy generated by some other particle 
instead of the heavy sterile neutrino. 



G. Summary of constraints 



III. MODELS WITH LOW SCALE TYPE I SEE-SAW 



Let us summarize this section. A theory where the DM ster- 
ile neutrino was in thermal equilibrium at some moment dur- 
ing the evolution of the universe should satisfy the following 
set of constraints: 

• From X/7-ray observations, we have the model inde- 
pendent upper limit on the radiative decay width of the 
DM neutrino N% (see Eq. (25)): 
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Note that this is a conservative value; c.f. footnote on 
page 5. 

This limit translates to the limit on the sterile-active 
neutrino mixing angle Eq. (24) and to the bound on 
the mixing between the left and right gauge bosons 
Eq. (28). 

• From the structure formation requirements (Ly-a 
bound), the mass of the sterile neutrino is constrained 
in the same way as the mass of a thermal relic, i.e. 



Mi > 1.6 keV . 



Let us start from the analysis of the models where the ac- 
tive neutrino masses are generated by a "type I" seesaw for- 
mula. This means that Ml = in the neutrino mass ma- 
trix (2). The mixing angles (7) are bounded from above by 
the requirements on the decay width of the sterile neutrinos — 
by the x-ray observations for the DM neutrino angle Q\ (see 
Eq. (24)), and by the long enough lifetime of the entropy gen- 
erating neutrino angle 82 (additional generation of the entropy 
by the third neutrino does not change conclusions). A conve- 
nient way to parameterise the Dirac mass matrix mo, sep- 
arating parameters in the active and sterile neutrino sectors, 
is provided by the Casas and Ibarra parametrization [35] re- 
viewed in Appendix C. Using Eq. (C9) for the Dirac masses, 4 
we get 



with 
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• The right abundance of the sterile neutrino can be then 
achieved by an out-of-equilibrium decay of a long-lived 
heavy particle. We will use another sterile neutrino of 
the model, N2, for this purpose, but most considerations 
here can be also applied to another long-lived particle 
present in the theory. 

To provide proper the entropy dilution, Eq. (14), N2 
should decouple while relativistic and has decay width 
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Here, R is a complex orthogonal matrix, describing the de- 
tails of the mixing between the sterile and active sectors, and 
it can be parameterized by three complex angles W12, 0J13, 
and cj 2 3 as in Eq. (C8). Let us check whether we can satisfy 
the bounds on the mixing angles if the active masses m, are 
consistent with the observed neutrino oscillation mass differ- 
ences, summarized below. The current best-fit and 3er ranges 
are (see Ref. [36]) 
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• At the same time, the heavy neutrino N 2 should de- 
cay before BBN, which bounds its lifetime to be shorter 
than approximately 0.1 -j- 2 s. Then, the proper entropy 
can be generated only if its mass is larger than 
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• The entropy is effectively generated by out-of- 
equilibrium decay (see Sec. II E), if the particle decou- 
pled while still relativistic. If this particle is one of the 
sterile neutrinos, then its decoupling happens at temper- 
ature (10), and it leads to the bound on the right-handed 
gauge boson mass 



In the following discussion we will for convenience order the 
active neutrino masses as m\ < r«2 < TO3. From Eq. (35), 
we get for the first two sterile neutrinos 



M\6\ = m 3 |sinwi3| +m 2 |coso;i3| |sinwi 2 | 



mi I cos wi3| 2 |coscji2| 2 



(38a) 



M 2 Q\ = m 3 |coswi 3 | |sino;23| 



+ m2|cosu;23 cos CJ12 — sincj23 sina;i 3 sinwi2| 

+ TOi|cosa>2 3 sinwi2 + sinw23 sinwi3 cosa;i2| 2 ■ 

(38b) 



M > 



1 

~T/a 
5*f 



M 2 \ 
lGevJ 



3/4 



(10 16) TeV 



(34) 4 As far as we are using in this section the basis with diagonal Mr, the 
right-handed mixing matrix is trivial, Vr = /. 
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Note that as far as we ordered the active neutrino masses, if 
we change mi to zero, and replace 7713 by m%, the right-hand 
sides of Eqs. in (38) can only become smaller. We can also 
confine ourselves to the real values of the mixing angles, as far 
as the sine and cosine absolute values only become larger for 
complex angles, and the inequality \z — w\ > ||z| — |io|| is 
used to transform the square of the difference of the angles in 
Eq. (38b). Thus, the following inequalities should be satisfied: 



Mi6\ > m 2 {sin 2 cj 13 + cos 2 w 13 sin 2 u; 12 } , 



(39a) 



M 2 9 2 > 77x2! cos 2 W13 sin 2 lu 23 + (|cos w 2 3||coso;i2| 

— |sinw23||sincji3||sintJi2|) 2 | . (39b) 

The minimum of the sum of the right-hand sides is 777,2, and 
therefore the following very simple inequality always holds: 



Mld\ + M 2 d\ > 777 2 > A777 so l 



(40) 



The second inequality is trivially fulfilled, since in all possi- 
ble mass hierarchies the mass of the second (in mass) active 
neutrino is larger than Am so i. The meaning of the inequal- 
ity (40) is very simple — one can not generate active neutrino 
masses with type I seesaw formula without sufficient mixings 
between the active and sterile neutrino sectors. Note in pass- 
ing that the cancellation is possible in another direction — one 
can have very small active neutrino masses and large active- 
sterile mixings. 

Now, we are ready to compare the requirement from the 
observed active neutrino masses, Eq. (40), and the DM bounds 
on the mixings. The angle 9 2 can be bound from the width 
required to generate sufficient entropy, Eq. (20). Estimating 
the width of the heavy neutrino as (see Appendix A) 



> 



G|Mf 

1927T 3 



we have 



M 2 0i < 1. 



10- 3 gi /2 



vy 



keVV 



(41) 



(42) 



It can be clearly seen that for all possible masses M\ and M 2 , 
this is much smaller than Am so i. 

The contribution of the DM sterile neutrino itself can be 
larger. From Eq. (24), we have 



Mi6\ < 1.8 x 10' 



lkeVV 



(43) 



Together with the Ly-a bound on the WDM mass, Eq. (16), 
this contribution again violates Eq. (40). Thus, we conclude 
that the small mixing angles, required by the proper DM abun- 
dance and good DM properties in the model, prevent genera- 
tion of the observed active neutrino masses by the type I see- 
saw formula. 5 



Note, however, that without the hy-a bound it would have been possible 
for very light WDM, with Mi < 1.2 keV. 



IV. TYPE II SEESAW— WORKING EXAMPLE 

In the previous section, we have seen that if one of the 
not DM-like sterile neutrinos is responsible for entropy pro- 
duction, it is impossible to get the observed active neutrino 
masses with the type I like seesaw. Here, we will present a 
working example of the sterile neutrino DM in the framework 
of a LR symmetric model, where the active neutrino masses 
are generated by the contribution of the type II seesaw. 

We will continue to work in the framework of the SU (3) x 
SU(2) L x SU(2) R x U(1)b-l model, sketched in Sec. II A. 
Here, we will concentrate on a properly LR symmetric model, 
where the left- and right-handed leptons are treated symmet- 
rically. One has the usual SM doublets ip % L ; i = 1, 2, 3, and 
in addition 3 right-handed neutrinos which form together with 
the 3 charged right-handed leptons the SU(2)r doublets i/j r . 
The Higgs sector consists of one SU(2) l triplet, one SU (2) R 
triplet, and one bidoublet. In such a model, the mass matrix 
for the neutrinos has the pattern 



M = 



Jlv l yv 
y T v f R v R 





f Ml 


m D \ 


)-' 


\ m D 


Mr ) 



(44) 



where the Majorana blocks on the diagonal come from the 

coupling of ijf L Cip J L and ip R Cip R with the triplets Al,r, 
respectively, and the Dirac-type ones from the coupling of 
^ipR with the bidoublet <f> and its complex conjugate <f> = 
T 2 <j)*T 2 . The VEVs of the neutral components in A^^ are 

called vl,r\ whereas, the SM scale v = \J k\ + K 2 = 
174 GeV is a combination of the bidoublet VEVs Ki and k 2 . 
These VEVs are related by the expression 



VLVR 



(45) 



where x is a function of the parameters in the Higgs poten- 
tial, which is naturally of order one (for more details, see e.g. 
Ref. [11]). 

In the following, we postulate exact discrete LR symmetry. 
In general, it can be realized in two different ways: as a C 
conjugation or as a parity symmetry. In the former case, it is 
required that ,/x = fn, V = y T , and this is what we use. 

With such a model, it is now possible to satisfy all the re- 
quirements from Sec. II. Let us consider the type II seesaw 
formula following from block diagonalization of Eq. (44) and 
the assumption 0(Mr) > <D(m D ) > 0(M L ): 



v 

VR 



vl/l - —yf R 1 y T 



(46) 



After applying the conditions of discrete left-right symmetry, 
y = y T and f = f L = f R , one arrives at 



Vhf 



vr 



-yf 1 y 



(47) 



To simplify the calculations, we further assume for illustration 
that the Dirac- Yukawa y is proportional to the triplet Yukawa 
/, i.e. y = pf, where p is a number. Equation (47) then goes 
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into 



m„ 



v 2 p 2 
vl / 

VR 



(48) 



In this case, all Yukawas are diagonalized by the same 
transformation — the transformation which brings m„ into di- 
agonal form, i.e. the PMNS matrix. The ratios of the eigen- 
values of the matrices on both sides of the equality are then 
the same 



mi 
m 2 



h 



Mi 

A/2 



(49) 



Thus, the mass spectrum of the sterile neutrinos (or, specifi- 
cally, the BBN requirement (22)) leads to the same hierarchi- 
cal active neutrino spectrum 



mi 
m 2 



< 5.9 x 10~ 7 



(50) 



This implies that the lightest active neutrino should be very 
light, and we can have either normal or inverse hierarchy. For 
definiteness, we will use the normal hierarchy for our exam- 
ple, though the inverse one works equally well (one should 
only take into account that in the latter case the M 2 ~ M3, 
r 2 ~ Ts and both Ni and 7V 2 generate the same amount of 
entropy). As far as the active neutrino mass hierarchy is fixed, 
we have m 2 ~ Am so i and m.3 ~ Am atm . We can then get 
the mass for the third sterile neutrino from 



Mo 



m 3 
m 2 



Mo 



(51) 



The active-sterile mixing angles in the case of proportional 
Yukawa constants are all the same and equal to 



(52) 



while the mixing angles for individual flavors are proportional 
to the PMNS matrix 



%i — [y al — • 

VR 



(53) 



Thus the decay width T 2 is proportional to 8 2 (see Appendix 
A). The value of 2 is then defined from the requirement of 
the sufficient entropy production, Eq. (20), and depends only 
on Mi and M 2 . 

At this moment the only free parameter left is the VEV ratio 
x, and everything can be expressed via x, Mi, M 2 , and m 2 ~ 
Am so i, 771 3 ~ Amatm. From Eqs. (48) and (45), we get 



vr 



(54) 



The VEV of the left-handed triplet A l is then given by 



Together with Eq. (54), Eq. (52) determines the proportional- 
ity constant p: 



P = ^Wr/v 2 . 
The full mass matrix (44) is then given by 



M = 



where 



U* 

u* 



M 



diag 
L 
diag 



M 



diag 
diag 



W 
W 



in 



M 



diag 
D 

diag 



^ M diag 

V L L 
VL 



l r 



i'r 



M^ ag 



and 



diag (Mi, M 2 ,M 3 ) 



(56) 
(57) 

(58a) 
(58b) 

(59) 



Let us fix now the input values. It was mentioned before that 
x ~ 0(1) is natural in the LR symmetric model, therefore 
we simply choose x = 1. For the masses of the DM and 
the entropy producing sterile neutrinos, we take the smallest 
possible ones (see Sec. II): 



Mi = 1.6 keV , 
Mi = 2.7 GeV . 

With this input, we obtain 

5.2 x 10" 9 eV 



(60a) 
(60b) 



mi 



m 2 



'sol 



= 8.7 x \{T i eV 



AmL = 4.9 x 10- 2 eV , 



-15 



m 3 c 

M 3 = 15.1 GeV , 
= 0| = 2.3 x 10 
v R = 9.67 x 10 4 TeV 
v L = 313 keV , 
p = 0.027 . 



(61) 



We also plot the values of Q\ and vr for several Mi and M 2 > 
2.4 GeV in the Figs. 3 and 4. Because of the smallness of 
2 compared to m 2 /M 2 and its suppression with M 2 or M 2 
(see Eq. (42)), vr given by Eq. (54) is effectively independent 
of 2 . Therefore, the curve of vr has only a very weak Mi 
dependence. However, for bigger Mi one has to account for 
the BBN bound (22) on M 2 . 

One can check that none of the bounds, summarized in 
Sec. II G is violated. Also the mixing angle Q\ correspond- 
ing to our DM neutrino is much lower than its upper bound, 
Eq. (24). However, we should also choose the Higgs potential 
to have very small mixing between the left and right gauge 
bosons (see Eqs. (28) and (29)). 

The right-handed scale vr is large, and the additional gauge 
and Higgs bosons are not observable (they all have masses 
oc vr). The famous p parameter 
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M 2 , GeV 



FIG. 3: 82 as a function of M.%. Mi = 1.6 keV (continuous); 
Mi = 5 keV (dotted); Mi = 20 keV (dashed); Mi = 40 keV 
(dashed-dotted). It is accounted for the lower bound, Eq. (22). The 
long-dashed (red) line shows the ratio (Mw/M) 4 and illustrates that 
processes mediated by Wr bosons can be neglected in the decay rate 
of N 2 . 





5 ; 


TeV 




m 
C 


3 ; 




2 ■ 
1 ; 



20 40 60 80 

M 2 , GeV 



DM sterile neutrino), (very) large number of degrees of free- 
dom in the early universe at the DM neutrino freeze-out, or 
subsequent dilution of its density by the out-of-equilibrium 
decay of a massive particle (another sterile neutrino). We fur- 
ther analyze this last possibility as being the most natural 6 
and formulate a set of requirements for this scenario. In short, 
these requirements bound the mass of the DM sterile neutrino 
from below from structure formation considerations, limit its 
mixing angle with active neutrinos and constrain mixing be- 
tween the SM (left) and additional (right) gauge bosons from 
the radiative decay of the DM sterile neutrino, fix the lifetime 
of the heavier sterile neutrino from the requirement of the di- 
lution of the DM abundance down to the observed value, and 
finally constrain the mass of these heavier sterile neutrinos 
from the big bang nucleosynthesis considerations. 

We demonstrated in this scenario that the type I low scale 
seesaw mechanism of generating masses for the active neu- 
trino can not lead to sufficient dilution of the DM abundance. 
At the same time, we provide a working example, where the 
active neutrinos are generated by a type II style seesaw in the 
context of an exactly LR symmetric theory. The provided gen- 
eral constraints and observations can serve as a basis for the 
search of a grand unified theory with WDM sterile neutrinos. 
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FIG. 4: vr as a function of Mi. The dependence on Mi is very 
weak. 



Appendix A: Decay widths of a sterile neutrino 



which is equal to 1 at tree level in SM also gets a negligible 
correction which is equal to 



2KI 



■4|« £ | 



(63) 



in the LR symmetric model [37]. For the small vl of the order 
of MeV, the deviations are well below the current experimen- 
tally allowed deviation of the order O(10~ 4 ) [38]. 



V. CONCLUSIONS 

In this paper, we analyzed the possibility to have a keV 
scale sterile neutrino warm dark matter in gauge extensions 
of the standard model. We found that it is possible to cir- 
cumvent the naive expectation of significant overproduction 
of dark matter in case of a light particle (sterile neutrino) de- 
coupling from the thermal equilibrium while still relativistic. 
The possible ways out include a low reheating temperature (so 
that the thermal equilibrium is never reached by the would be 



In the mass range 2.7 GeV < Mi < Myy, a sterile neutrino 
N2 dominantly decays into leptons and spectator quarks. The 
corresponding partial widths can be calculated in the vMSM, 
because additional boson interactions, which usually appear 
in more complicated models, are of high scale — of order 
0(M) — compared to that of the electroweak scale. To make 
use of the i^MSM results, we have to compare, for being on 
the safe side, the suppression factors in the z/MSM ~ 0{\ 0i | ) 
with that of additional interactions ~ 0(Mw /M) 4 which ap- 
pear in the models we are interested in. Furthermore, the mix- 
ing of the new bosons should be small compared to 
otherwise, there could be significant contributions from pro- 
cesses where new bosons mix with the SM ones (see, for ex- 
ample, Appendix B). These effects are neglected in the fol- 
lowing calculations. One can see that for most practical pur- 
poses it is the case, as far as the bound on the gauge boson 



6 Another natural possibility is achieved in the i^MSM model [6, 7], where 
the sterile neutrinos are the only extension of the SM, and then the keV 
sterile neutrino does not enter thermal equilibrium up to Planck scale tem- 
peratures. 
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mixings (28) is much stronger then those for the active-sterile 
neutrino mixings. 

Moreover, these additional contributions do not affect the 
conclusions in the main part of the article. Really, in Sec. Ill 
additional interactions can only result in a stronger bound and 
therefore the conclusion remains the same. In the type II see- 
saw model discussed in Sec. IV, we need the exact value of 
the width. However, if we assume no mixing of the W bosons 



(C = 0), the contributions of the Wr boson mediated pro- 
cesses are in the considered mass range negligible small (cf. 
Figure 3). 

In the following, we give all relevant formulas for the decay 
rates (at tree level) of a sterile neutrino N 2 with a mass M 2 
above the BBN bound 2.7 GeV (cf. (22)) and below the SM 
W-boson mass M w ~ 80 GeV [39]: 



with 



r 2 [n 2 -> i-^ffivfi 
r 3 (n 2 -> vjpp 



G 2 F Mj 

1927T 3 



1927T 3 

G 2 F Mj 

1927T 3 



£|<U 2 , 



■ \02 a \* {l - 8xf + 8xf -xf -Uxflogxf) , XI 



(Ala) 



max [Mi a , Mi p ] 



M 2 



(Alb) 



(Ci • (1 - 6 a p) + C 3 ■ 8 aP ) {{l - Ux 2 - 2xf - 12s?) yf\ - Ax 2 



+ 12xf (xf - 1) L + 4 (C 2 • (1 - M + C* 4 • 6 a p) [xf(2 + Wxf - Uxf) J 1 - 4^ 



+ 6xf (1 - 2xf + 2xf) L 



(Ale) 



and 



L = log 



3xf - (1 - x 2 ) ^Jl-Ax 



x 2 (1 + vT=I£f) 



d = - (1 - 4 sin 2 6*^+8 sin 4 9 W ) , 



C 3 = - (1 +4sin 2 



xi 



M 2 



C 2 = -sm 2 9 w (2 sin 2 6 W - l) , 
C 4 = - sin 2 9 W (2 sin 2 6 W + 1) . 



The formulas for the decay modes into quarks are presented 
below. In the range 2.7 GeV < M 2 < Mw, it is sufficient 
to use the free quark approximation for the decay products. 
We give these formulas in the approximation where M 2 is 
much heavier than the decay product masses (unlike above 
(Eqs. (Al)) for the lepton decays). The corrections are impor- 
tant at the threshold, when new decay channels open. How- 
ever, at high mass M 2 , this introduces a rather small relative 
error, because the number of open channels into light particles 
is significant and provides the main part of the decay width. 
The exact analysis would smooth the discontinuities of the de- 
cay width at the mass thresholds (Fig. 5): 

T 4 (JV a l-UD) - • 3 • \V UD \ 2 ■ \0 2a \ 2 , (A2a) 

C 2 M 5 

T 5 (N 2 -> v a qq) = -^-f ■ 3 • S« • \6 2a \ 2 , (A2b) 



with 

z q = (al) 2 ■ ((si) 2 + {g q R ) 2 ) ■ 

The factor 3 is the color factor and Vjjo are the Cabibbo- 
Kobayashi-Maskawa-matrix elements. The coupling con- 
stants cjl and gn correspond to the coupling of the Z boson 
to left- or right-handed particles, respectively. For a fermion 
/ with weak isospin component l( and charge qf, one has 

9l = 4 -If sin2 8 w 1 (A3 a) 

9r = ~Qf sin 2 9w ■ (A3b) 

Table I gives the required values of the charges. In the mass 
range of M 2 mentioned above, the Majorana neutrino total 
decay rate T^ 2 is a sum of all rates presented above multiplied 
by a factor of 2, which accounts for charge-conjugated decay 
modes. 
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fermions 


9l 


9r 


U = u,c,t 
D = d,s,b 


9l = \ 

U _ 1 2 2 
9h — 2 3 S W 

f) D - 1 -1- 1 s 2 

9h — 2 ' 3 S W 


9r = 

n U - 2 e 2 
9r — 3 S W 

n D - 1 <i 2 



TABLE I: Coupling constants 



When M 2 exceeds the SM Z-boson mass Mz =91 GeV 
and if contributions of new interactions are negligible, the 
sterile neutrino predominantly decays into a SM gauge boson 
and a lepton. Then its total decay width is given by 



N 2 



G F Mj 
1 

+ 2 



2M W 
Mf 



2M| 
~W 



1 - 



Mf 
Ml 



Ml 



E 




(A4) 



0.01 0.1 

M 2 , GeV 



FIG. 5: The ratio rjv 2 /2ri calculated in the specific LR model of 
Sec. IV. We used Eq. (Al) together with the formulas for the two- 
body decays into mesons [39] for the left (red) curve and Eq. (Al) 
together with Eq. (A2) for the right (blue) curve. The dotted curves 
correspond to the region, where both approximations are not entirely 
reliable. 



In between, where Mw < M2 < Mz, one can approximate 
the width by the first term in Eq. (A4). 

Below M2 ~ 2 GeV, it is important to consider mesons 
instead of quarks as final states for the sterile neutrino de- 
cay. Therefore, instead of the three-body decay modes into 
spectator quarks (A2), one has to use the corresponding two- 
body ones into mesons; see [39] for the decay width formulas. 
In our case, this mass range of M2 is forbidden by the BBN 
bound (22) and therefore we do not list them here. Neverthe- 
less, to get a feeling of the behavior of the total width Fjv 2 , 
we show in Fig. 5 the ratio rjv 2 /2ri calculated in the specific 
model described in Sec. IV, using both the free quark and chi- 
ral meson approximations. It is clearly seen that the transition 
between two approximations happens around 1 GeV. Because 
of Eq. (53), the total decay width is proportional to B\ and 
therefore the plotted ratio is independent of this quantity. In 
more general models, Eq. (53) will no longer be valid. How- 
ever, as one recognizes by considering the formulas together 
with the definition of O2 (cf. Equation (7)), there will be no 
significant difference, especially for heavy masses M2, so that 
Fig. 5 can be used as a good estimate in such models. Note 
that in the region M 2 ~ 0(1) GeV (dotted in Fig. 5), de- 
cays into spectator quarks more and more replace decays into 
mesons and therefore one has to carefully reanalyze the given 
formulas, if one is interested in this mass range. 



vector of the outgoing photon. The Ward identity for the elec- 
tromagnetic current implies that q^M.^ must be zero; there- 
fore Ar M must have the form 



[ia^q y (a L L + a R R)] Ui , 



(Bl) 



(1 - 75 )/2 and R = (1 



where a^ v = {i/2)[y», 7"], I 
75)/2 are the projectors of chirality. <tl and ctr are numerical 
coefficients with dimension of inverse mass. The partial decay 
width for N\ — > 1^7 is then given by 



[Ml - mf) 3 
IQirMl 



(tl\ +Wr 



(B2) 



By comparing the Lagrange term for the charged current (1) 
combined with Eq. (8) and the transformation rule which di- 
agonalizes the neutrino mass matrix in Eq. (2) 



Appendix B: Radiative decay width 

Here, we give some details of calculation of the width for 
the radiative decay Ni — > shown in Fig. 2. We will follow 
Ref. [33], where general formulas for this type of process are 
given. In our case, N\ denotes a heavy sterile neutrino with 
mass Mi, Vi one of the active neutrinos with mass and 7 
a photon. The neutrinos are considered as Majorana particles. 

The amplitude for such a decay is ee*^(q)Ai 11 , where e is 
the electric charge of the positron and e*(q) the polarization 



fa) = ( A *)(»") , (B3) 

with that given in Chapter 5 of [33], we can easily calculate 
the coefficients and ctr. Supposing from the very begin- 
ning, that the right-handed scale is much larger than the left 
one, M ^> Mw — 80.4 GeV, and neglecting the active neu- 
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trino masses, we get 



2 

9 e 



32 JV/2, 7T 2 



^ {cosCsinCA^D^m^J-^) 
+ cos 2 C^B al A/ 1 F(r Q )} , 



(B4a) 



2 

3 e 



"32Af2,7T 2 



^ { cos C sin C C ai B„i m, a J(r fl )} , (B4b) 

a—e,fi,T 

where F{r a ) and .F(r ) are functions of r a = m 2 /Myy- In 
our case, we have in good approximation F(r a ) ~ —3/2 and 
^(ra) ~ 4. The exact expressions for these functions were 
calculated by us and do agree with that given in Ref. [40]. 

Because of the Majorana nature of our ingoing and out- 
going neutrinos, we also have to add the contribution of the 
complex conjugated process to our amplitude. This is easily 
obtained out of Eq. (B4) by putting in the substitutions 



A,B^A*,B* and C,D -> C*,D* , 
75 -> -75 => L, R -> R, L , 



(B5) 



and an overall negative sign coming from the photon vertex. 
After adding the derived crj, and <tr, it is easy to see that 
Wl\ 2 = Wb,\ 2 , where 



9 

g e 



32 Myf ir 2 



a=e,jU,r 



4 cos C sin C ^2 (A al D al - C ai B al )mi a 



3 2 

- cos 
2 



(B6) 



By putting this into Eq. (B2), we obtain 

G 2 F aMl 



Ni—tjVi 



64 tH 



4 cos CsinC ( A aiD a i-C ai B al ) 



-cos 



\a — e.f-i.r / 



(B7) 



7 Note that our results do not coincide with the formulas in [40]. This is be- 
cause of a mistake in the second term of the third line of Ew. (10) in [40]. 
The correct labeling of the transformation matrices should be P a BQaA 
instead of P a AQaB- In our notations, where a sterile neutrino (with mass 
eigenstate index 1 ) decays through the radiative process into an active neu- 
trino (with mass eigenstate index i), the expression P a BQaA translates 
into BaiC a i which is contained in Eq. (B4b). 



Here, Gf is the Fermi constant, a is the fine-structure con- 
stant, and mi a is the mass of the charged lepton propagating 
in the loop. 

The total width of the radiative decay is given by 



(B8) 



In a model where a seesaw mechanism of type I or type II is 
responsible for the small active neutrino masses, the transfor- 
mation (B3) is given by Eq. (3). Putting this into our formulas, 
we get out of Eq. (B7) the expression (26). 



Appendix C: Casas-Ibarra Parametrization 

In this part of the Appendix, we describe the approach of 
parametrizing the Dirac- Yukawa matrix, which was proposed 
by Casas and Ibarra [35]. Here, we want to give a short review 
of the generalised version which also applies to the type II 
seesaw mechanism [41]. 

Let us consider a Majorana mass matrix with the pattern 



( M L 


m D 




V m D 


Mr ) 





Jlv l yv 
y T v f R v R 



(Cl) 



The type II seesaw formula can be written in the form 

m v - Ml — —mDM^m^) , (C2) 

where m v is the active neutrino mass matrix. Let us define the 
symmetric and in general complex 3x3 matrix 



X,, = 771 1. 



M 



L ■ 



(C3) 



This matrix and Mr can be diagonalized by unitary transfor- 
mations: 



x v = v;x^v} = \v;{xf a ^ 

Mr = V^M^Vt . 



R 



(C4a) 
(C4b) 



Multiplying Eq. (C2) by 

by l\v:(x^s)f 7 



V*(X^ 



from the left and 



from the right and using Eq. (C4), 



we find 



with 



I = RR J , 



(C5) 



R = ±i(X^s) -*Vjm D V R (i< ias ) *. (C6) 

Equations (C5) and (C6) mean that the type II seesaw relation 
requires R to be a complex orthogonal matrix, but otherwise 
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does not constrain it. In this way we obtain for the Dirac-type 
Yukawa coupling in the basis where Mr is diagonal 

m D = vy = ±i V: \fx<F s R yj ' M^ S V^ , (C7) 

where R is an arbitrary complex orthogonal matrix. It can be 
parametrized as 

R = ± i?i2 i?i3 R23 , (C8) 

where Rij is the matrix of rotation by a complex angle u>ij in 
the ij plane. This is the so-called Casas-Ibarra parametriza- 
tion of the Dirac- Yukawa [35]. Note that this parametrization 
has its origin in the difference of the number of high energy 



and low energy parameters. There are less low energy param- 
eters, because the high energy ones are integrated out. The 
latter cannot influence the low energy theory, and therefore 
can be parametrized arbitrarily. 

The formula for the type I seesaw can easily derived out of 
(C7). Because of Ml = 0, X v corresponds in this case to 
the active neutrino mass matrix m v . Thus the basis transfor- 
mation matrix V v in Eq. (C4a) is the PMNS matrix U and we 
arrive at 



m D = vy = ±i U*\] mf ag R \J M^ &e V^ . (C9) 
In chapter III we make use of this parametrization. 
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